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We study a tight binding model including both on site disorder and coupling of the electrons to ran- 
domly oriented magnetic moments. The transport properties are calculated via the Kubo-Greenwood 
scheme, using the exact eigenstates of the disordered system and large system size extrapolation 
of the low frequency optical conductivity. We first benchmark our method in the model with only 
structural disorder and then use it to map out the transport regimes and metal- insulator transitions 
in problems involving (i) scattering from random magnetic moments, and (tt) the combined effect 
of structural disorder and magnetic scattering. We completely map out the dependence of the d.c 
conductivity on electron density (n) the structural disorder (A) and the magnetic coupling (J'), 
and locate the insulator-metal phase boundary in the space of n — A — J'. These results serve as a 
reference for understanding transport in systems ranging from magnetic semiconductors to double 
exchange 'colossal magnetoresistance' systems. A brief version of this study appears in our earlier 
paper Europhys. Lett. 65, 75 (2004). 



I. INTRODUCTION 

The most commonly studied case of localisation per- 
tains to non interacting electrons in the background of 
structural disorder. There is a large body of work 1-4 , 
analytical and numerical, as well as experimental stud- 
ies, that have focused on this problem. The principal 
qualitative result of these investigations is that in one 
and two dimensions all electronic eigenstates are localised 
for arbitrarily weak disorder, while in three dimension we 
need a critical disorder for complete localisation. In three 
dimension, at a given disorder, all states beyond an en- 
ergy e c of the band center are localised and the system 
is metallic or insulating depending on whether the Fermi 
level, cf , lies in the region of extended states or localised 
states. The 'mobility edge', e c , collapses to the band 
center as the disorder is increased, driving the Anderson 
metal- insulator transition (MIT). 

The presence of magnetic moments in a metal brings in 
several new effects, depending on the strength of electron- 
spin coupling (J'), the concentration of moments (n mag ), 
the extent of disorder, and the 'character' (small or large 
S) of the moment. 

In the 'quantum limit', 2S *~ 1, and for antiferro- 
magnetic coupling, the basic physics is contained in the 
Kondo effect. For n mag <C 1, the magnetic moments act 
as 'Kondo impurities' whose effects 5 are now well un- 
derstood. For n mag ~ 1, i.e, the concentrated Kondo 
limit, there can be various phases depending on electron- 
spin coupling and disorder. The ground state could be 
a non magnetic 'heavy Fermi liquid' 6 , or a spin glass 7 , 
or a magnetically ordered state 8 . The physics of these 
Kondo lattice, with quantum spins, is a vast area of re- 
search. In this paper, however, we will avoid the issues 
of heavy fermion physics and focus instead on electron- 
spin systems involving "large 5", i.e, effectively 'classical' 



moments. 

For classical moments also, the effects vary depending 
on n magi electron density, J', and the extent of disorder. 
A wide variety of magnetic systems 9-17 are described, to 
a first approximation, by electrons locally coupled to d or 
/ moments, with 2S 3> 1, and moving in a structurally 
disordered background. The magnetic ground state could 
be ferromagnetic, or a more complicated ordered state, 
or a spin glass. Transport often involves insulator-metal 
transitions and colossal magnetoresistance. The simplest 
Hamiltonian capturing these effects is: 



H = -t 



C icr C i<* 



s„ (1) 



The t are nearest neighbour hopping on a simple cubic 
lattice. The random on site potential, e^, is uniformly dis- 
tributed between ±A/2. The sites R„ are a subset of the 
cubic lattice sites, R$, and correspond to the magnetic 
'dopant' locations. Even with this simple model there are 
four dimensionless parameters in the problem. These are 
disorder A/t, magnetic coupling J'S/t, electron density 
n (controlled by /z), and the 'density' of moments n mag . 
We will eventually study the n mag — I case, but retain 
a more general structure right now. We absorb S in our 
magnetic coupling J 7 , assuming |S*| = 1. 

Real materials have band degeneracy and additional 
interactions but the basic physics of several currently in- 
teresting materials arise as limiting cases of the model 
above, (i) The II- VI diluted magnetic semiconductor 9-11 
(DMS) Gai^Mn^As, exhibiting high ferromagnetic T c , 
correspond to n mag <C 1, J'/t ~ I, weak disorder, and 
low electron density, n < n mag . (ii) The Eu based mag- 
netic semiconductors 12,13 , EuBg etc, involve n mag = 1, 
since every Eu atom has a moment, J'/t 1, low carrier 
density, and possibly weak disorder. (Hi) The 'colos- 
sal magnetoresistance' (CMR) manganites 14 , specifically 
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Lai-^Sr^MnOs, involve n mag = 1, J'/t 3> 1, high elec- 
tron density, and moderate 'effective disorder'. To de- 
scribe the more strongly resistive manganites, the Ca 
doped systems, say, one requires additional electron- 
phonon interactions. (iv) The amorphous magnetic 
semiconductor 15,16 , a-Gd^Sii-u;, corresponds to J'/t 3> 
1, A/t > 1, and n mag ~ n ~ 0(0. 1). Finally, (v) the 
traditional metallic / electron magnets 17 , correspond to 
n-mag = 1, and moderate to strong J' . 

The focus in the materials above is often on magnetism 
rather than localisation effects. However, many of them 
have rather large resistivity in the paramagnetic phase, 
and a-GdSi, for example, shows a metal-insulator tran- 
sition at T = itself, on lowering carrier density. Since 
there is no direct spin-spin interaction in these systems, 
the local electron-spin coupling controls both the magnetic 
properties and the character of the electronic state. 

The intimate coupling between charge transport, local- 
isation effects, and magnetism in these systems suggest 
that we need to look beyond the traditional boundaries 
separating 'magnetism' from transport and localisation 
studies. A complete study of electronic resistivity as a 
function of temperature, for any of the materials above, 
requires a solution of the magnetic problem first. Since 
the moments are assumed to be classical, the electrons 
can be imagined to move in a static background compris- 
ing the (quenched) structural disorder and annealed spin 
disorder. Evaluating the distribution of the annealed dis- 
order is a non trivial problem, particularly in the strong 
coupling (large J') context that is experimentally rele- 
vant. We will touch upon this in the next section, but 
this paper is concerned with transport and localisation 
effects in the fully spin disordered phase. In this limit, 
we will present a comprehensive discussion of the resis- 
tivity arising from the interplay of structural disorder 
and 'paramagnetic' scattering, and map out the metal- 
insulator phase diagram in terms of electron density, dis- 
order and magnetic coupling. 

There have been some studies of electronic transport 
in the background of random spins and structural dis- 
order, acting independently or together. Among these, 
the Anderson localisation problem itself has been exten- 
sively studied, via perturbation theory 18 , self-consistent 
schemes 19 , numerical techniques 20 , and mapping to a 
field theory 21 . Most of the qualitative issues in this con- 
text are essentially settled. Weak magnetic scattering in 
a structurally disordered system has been studied 22 ' 23 in 
the early days of weak localisation (WL) theory to clar- 
ify the 'dephasing' effect of electron spin flip on quan- 
tum interference. In the opposite limit of strong cou- 
pling, corresponding to double exchange, localisation ef- 
fects have been studied 24 considering both magnetic and 
structural disorder. 

These efforts still leave a large and interesting part of 
A — J' space unexplored. To give a few examples, there 
is no discussion of the following: (i) the resistivity from 
purely magnetic scattering, as J' rises through the per- 
turbative regime to double exchange: this is the classic 



problem of paramagnetic scattering in 'clean' magnets, 
studied earlier at weak coupling 25 ' 26 , (ii) the effect of 
spin disorder on the Anderson transition, i.e, how the 
'anti-localising' effect of spin flip scattering, at weak dis- 
order, evolves into an insulator-metal transition (IMT). 
This is an instance of Anderson transition with bro- 
ken time reversal symmetry, and (Hi) the wide 'middle', 
where the effect of neither A nor J' is perturbative and 
their contribution to the resistivity is not additive (i.e, 
violates Mathiessens rule) . This is the regime relevant to 
DMS, CMR materials, and amorphous magnetic semi- 
conductors. 

The next section describes the transport calculation 
in detail. Following that we present results on transport, 
successively, in the structural disorder problem, the mag- 
netic disorder problem, and the simultaneous effect of 
both. This paper follows up on our earlier short paper 27 . 

II. COMPUTATIONAL SCHEME 

Although we will work with random uncorrelated 
spins, viewing the magnetic disorder as quenched, let us 
highlight how the 'true' spin distribution can be evalu- 
ated, and the limit where the background can be consid- 
ered random. Following that we describe our transport 
calculation method. 



A. The spin distribution 

The 'structural' variables ej are quenched, and have 
a specified distribution, but the spin orientations Si are 
not known a priori. The system chooses a spin configu- 
ration, at T = 0, to optimise the total energy. To calcu- 
late the 'true' ground state properties, or finite tempera- 
ture transport, we need to solve for the spin distribution 
first and then evaluate electronic properties in these spin 
background. Denoting the full spin configuration as {S^}, 
the spin distribution -P{Si} is given by: 

P{Sj} = Z- x Tre- (m 
Z = J VSiTre'l 311 

where Z is the full partition function of the system, 
and the 'trace' is over fermionic variables. Equivalently, 
the effective classical 'Hamiltonian' controlling the Boltz- 
mann weight for spins is: 

Heff{Si} = —log Tre-? H (2) 

H e ff is the fermion (free) energy in the background {Si}. 

To make more sense of the formal expression above, 
consider J'/t <C 1. In this case we can expand the 
fermion (free) energy in powers of J'. For a non 
disordered system this leads to the standard RKKY 
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coupling 28 between the classical spins, while the presence 
of structural disorder, leads to a 'bond disordered' RKKY 
model: H e ff ~ ^ f . JySj.Sj, where the exchange are 

~ J' 2 Xij the being the non local spin response func- 
tion of the disordered, J' = 0, electron system. Having 
obtained the effective spin Hamiltonian, the transport 
properties are to be calculated by considering electron 
motion in the backgrounds {e^S^} where the {S^} are 
equilibrium configurations of H e f f for a specified realisa- 
tion of disorder {ei}. 

At strong coupling, i.e, large J', the fermion trace can- 
not be analytically evaluated, and it is no longer possible 
to write an explicit spin Hamiltonian. We need special 
techniques to anneal the spins. The magnetic order and 
the complete transport properties in such (disordered) 
Kondo lattice models is discussed elsewhere 29 ' 30 . 

The complications of the magnetic problem can be 
avoided if we assume a spin distribution. The simplest 
distribution one can assume corresponds to uncorrelated 
random spins. This is physically relevant in two limits. 

(i) At sufficiently high temperature, compared to the 
magnetic ordering scales in the problem, the spins are 
essentially randomly fluctuating, with only short range 
correlation. The magnetic ordering scale for J'/t <C 1 
is ~ fi(n)J' 2 /t, while for J'/t ^> 1 the ordering scale is 
~ f2(n)t, where /i and ji are electron density dependent 
dimensionless coefficients and f™ ax ~ 0.1. Compared to 
the typical Fermi energy, ~ zt, where z is the coordi- 
nation number of the lattice, these scales are all small. 
We use a T = formulation for transport, i.e, we do 
not use Fermi factors, but given the smallness of T c /cf, 
our results would be relevant even in the 'real' param- 
agnetic phase, (m) Another situation in which a random 
spin configuration is appropriate is a 'spin glass', likely 
to occur in strongly disordered systems 31 . There are al- 
ways short range correlations in a spin glass but if we 
ignore their effect on transport then at all temperature 
the transport in such a system can be described, approx- 
imately, in terms of a random spin background. 



B. Conductivity calculation 

In the linear response regime, the Kubo formula 
can be used to calculate the conductivity of a system. 
The general expression 32 , involving matrix elements be- 
tween many body states, simplifies significantly for non- 
interacting systems. This 'Kubo-Greenwood' result can 
be computed purely in terms of single particle cigenfunc- 
tions and energies. 

The numerical difficulty with this result lies in imple- 
menting it for a finite size system, where the spectrum is 
discrete, with gaps 0(1/N), with N being the number of 
sites in the system. Since the 'd.c' conductivity involves 
transitions between essentially degenerate states at , it 
cannot be calculated with control on small systems. As 
a result, instead of computing the 'Kubo conductivity' 



it is more usual to compute the 'Landauer conductance', 
G, of a finite system coupled to leads 33 , and infer the 
conductivity from the length dependence of G. 

We pursue the Kubo approach, indirectly, through a 
calculation of the low frequency optical conductivity for 
a Lt x Lt x L geometry. The advantage of calculating 
the conductivity this way is, (i) it ties in with diagonal- 
isation that one may have to do for the magnetic prob- 
lem, and (m) it allows access to the optical conductivity, 
without added effort, and can reveal the significantly non 
Drude nature of a(ui) at strong disorder. The principal 
disadvantage is, this scheme cannot be pushed beyond 
TV <~ 10 3 — 10 4 , and is therefore not useful for accessing 
critical properties. 

For disordered non interacting systems, the Kubo for- 
mula, at T = 0, is: 

^H^^EK-^t^^-^-^)) (3) 

with A = 7re 2 /hao, ao being the lattice spacing, and n a = 
0{pL — e a ). The / Qj g are matrix elements of the current 

operator j x = itJ2i.cr( c l+xa ,cr Cl ^ ~ h - c )' between exact 
single particle eigenstates \ip a ), \4>p), etc, and e a , ep are 
the corresponding eigenvalues. 

The conductivity above is prior to disorder averaging. 
Notice that the 6 function constraint cannot be satis- 
fied for arbitrary frequency in a finite system. So we 
can neither calculate the d.c conductivity, <Jd c , directly, 
nor estimate cr(ui) at some arbitrary externally specified 
frequency. However, we can still calculate the 'average' 
conductivity over a frequency interval Alj, defined below, 
and we use the following strategy to extract a de- 
li) The average of <j(lo) over the interval [0, Aoj] is de- 
fined as 

a av (AtJ, fi,N) = — / a(oj,n,N)du) (4) 
^ Jo 

Alo can be set independent of N, but we will relate them 
by fixing: Auj = B/N. We fix B by setting Auj = 0.04 for 
N = 1000. The mean finite size gap is 12/1000 ~ 0.01, 
in 3d, much smaller than Alo. 

(ii) a av is averaged over N r realisations of disorder, to 
generate a av (Aw, fi, L). The 'noise' in a av (Au), fi, L) falls 
slowly, as l/\fW^. We use N r <~ 100 for the largest sizes, 
to keep the computation reasonable, and use a filter to 
smooth the data over a small window in /i. 

(Hi) We study the a av (Au, fi, L) for Lt = 6 and the 
sequence L = 24 to L = 64 in increments of 8 and ex- 
trapolate, Tcaidp) = liTni / ^ 00 a av (Au>, n, L). As speci- 
fied before, Auj = B/N. 

To calculate the full, disorder averaged, optical con- 
ductivity we use the inversion: <j(oj) = a av (ui) + ui . 
The cr(u)) results in this paper are mostly based on a 
6 x 6 x 32 geometry. 
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III. TRANSPORT IN THE ANDERSON MODEL 



The metal-insulator phase boundary and the criti- 
cal properties near the transition have been extensively 
studied 18 ^ 21 in the Anderson model. However, the actual 
resistivity seems to have received much less attention. As 
recently pointed out by Nikolic and Allen 34 , there is a 
wide regime in A, between the Born-Boltzmann end and 
the scaling regime, where there are no analytic theories 
of transport. We study this 'old problem' in some detail 
because the wealth of existing results provides a bench- 
mark for our method. There are very few exact results 
with which we will be able to compare our data in the 
magnetic scattering problems. 



A. Global features 
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FIG. 1. Panel (a). Variation of conductivity with u, and 
panel (b). density of states, for several values of A. 



The 'global features' of transport and localisation in 
the Anderson model are contained in Figs. 1-3. The data 
is obtained via the extrapolation procedure described 
earlier. Fig.l.(a) highlights the suppression in conduc- 
tivity with increasing disorder, across the entire band. 
The 'weakest' disorder in this case, A = 4, is probably 
already outside the semiclassical Boltzmann regime. A 
naive A 2 scaling of the resistivity still works, at the band 
center, between A = 4 and A = 8, but the same extrapo- 
lated to A = 16, would imply p(16)/p(4) ~ 16, while the 
ratio is actually ~ 10 3 . This figure captures the expected 
crossover from moderate scattering, roughly following 
Boltzmann scaling, to localisation as A — > A c <~ 16.5, 
the critical disorder 35 at the band center. It also pro- 
vides a glimpse of how the 'mobility edge' moves with 
increasing disorder, better quantified in Fig. 3. Note that 
for data at a specified system size, L = 16, 32 etc, shown 
later, the notion of a 'mobility edge' does not make sense, 
and all we observe is a crossover from small to large con- 
ductivity as /j, is varied. The change in a(fM, L) with L, 
and the L — > oo extrapolation, is crucial for identifying 
the mobility edge. 

The DOS plot, Fig.l.(b), illustrates the band broad- 
ening, and in Fig.3.(b) we have compared our band edge 
energy with earlier CPA results 36 . The (algebraic) aver- 
age DOS is featureless and non critical and does not play 
an interesting role in the problem. 

Since the band broadens significantly with disorder, 
<t(/z) by itself does not provide the density dependence 
of the conductivity. Fig. 2 takes into account the shift in 
jj, required to maintain constant density (with increasing 
disorder) and shows a(n). Our density is defined as av- 
erage number of electrons per site, so n max = 2. Since 
the model is particle-hole symmetric we show only the 
regime n = [0, 1]. 




FIG. 2. Variation of conductivity with carrier density, for 
several A, constructed from the o{u) and N(e) data in Fig.l 



i i i | i i i | i i i | i i 

(a) 




FIG. 3. Main panel, (a), shows the variation in mobility 
edge with disorder. We compare our results, circles, with 
earlier work 36 , triangles. Inset (6) shows the 'band edge', and 
(c) the fraction of localised states at large A. 
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To get a feel for the magnitude of the conductivity, 
which we measure in units of ire 2 /(frao), note that the 
Mott 'minimum metallic conductivity', (JMott, at the 
band center 4 , is roughly <~ 0.03e 2 /(?iao). Our dimcn- 
sionless conductivity a ca i c , shown in the figures, can be 
converted to real units, a 'actual, by using 

factual ~ 100 * (JMott * Ocalc 

where we use OMott = 0.03e 2 / (ha ). The results we show 
in the present spin degenerate problem includes a factor 
of 2 to account for the two spin channels. This is im- 
portant to compare with the magnetic scattering prob- 
lems later. The conductivity per spin channel falls below 
w 10" 2 for A > 8. This implies that beyond A w 8, 
a < (JMott in the Anderson model. 

The main panel in Fig. 3 shows the variation in mobil- 
ity edge with increasing disorder. Our error estimates are 
based on the shift in e c as we change from moderate to 
zero filtering of the tr(/x, L) data. We show some earlier 
standard result 36 for comparison. The best current result 
on A c is 16.5, our method yields A c ~ 17. Our results on 
the band edge, Fig.3.(b), match reasonably with earlier 
CPA based results. Note that while the mobility edge has 
a 're-entrant' behaviour, the fraction of localised states 
in the band, Fig.3.(c), increases monotonically with dis- 
order. 



B. Transport regimes 



There are tentatively three transport regimes in the 
Anderson model. These are (i) the perturbative Born 
scattering regime, described by the Boltzmann transport 
equation and the low order corrections in (kpl) -1 ■ This 
corresponds to A/W -C 1, where W — 12t. (ii) The wide 
intermediate coupling regime A/W <~ 0(1), and (iii) the 
'scaling' region, A — ► A c , near the MIT. 




FIG. 4. Variation of the resistivity, normalised to the Born 
resistivity, with disorder. The firm lines are fits of the form 
1 + ai(n)A 2 + o?2(n)A 4 . Inset shows the coefficients a.\,a.i. 



We analyse the data, with increasing A, in the se- 
quence (i) — ► (ii) — > (iii). 



1. Perturbative regime 

To leading order, the scattering rate from the disorder 
is r^ 1 ~ 2ttN(€f)(^1)- The second moment of the 
random potential is, (e 2 ) ~ A 2 /12. Since N(e F ) ~ 0.13, 
at the band center, Fig.l.(b), the scattering rate, Ta = 
t^ 1 « A 2 /(15t). 

The three related quantities which define Boltzmann 
transport are (i) the scattering rate, Ta, defined above, 
(ii) the (inverse) mean free path ao/l ~ 0.03(A/t) 2 at the 
band center, and (iii) the Born-Boltzmann conductivity 
<tb ~ 1.62(7re 2 /hao) * (t/A) 2 per spin channel, again at 
the band center. In addition the optical conductivity 
should have the Drude form <j(uj) = <J B /(1 + (uj 2 /T\)), 
and the 'width' in the optical conductivity can be checked 
against the magnitude of d.c conductivity. 

Using the form for as, the conductivity at A = 4, 
assuming Boltzmann transport, would be approximately 
0.1(7re 2 /?iao) per spin channel, i.e. ~ 0.2 in our units in- 
cluding spin degeneracy. Our data, Fig. 2, gives a value 
~ 0.19 at the band center. The crude Boltzmann scal- 
ing is obvious from the moderate A results in Fig.l and 
Fig. 2. In Fig. 4 we attempt to quantify the corrections to 
the Boltzmann result, still staying far from the localisa- 
tion regime. 

The weak localisation corrections that arise beyond 
Boltzmann transport control the resistivity in one and 
two dimension. These are quantum interference effects, 
formally arising from the 'Cooperon' corrections. A sim- 
ilar argument would lead us to believe that in three 
dimension 1 the leading correction beyond the Boltzmann 
results should be 8a oc — (kpl)^ 1 . Since (kpl)^ 1 oc A 2 
and the Boltzmann conductivity <jb oc kpl, the net 
conductivity would be expected to have the form a ~ 
a B (\ - 0((k F l)~ 2 )), i.e, (j(A) ~ A- 2 (l - 0(A 4 )). In 
that case, the resistivity should have a form p(A) <~ 
PB (A)(1 + 0(A 4 )). 

Fig. 4 shows p(A)/pb(A) plotted against A 2 for three 
densities. We avoid too low a density to keep the scales 
comparable. The data are fitted to p/pb = l+ai(n)A 2 + 
<22 (n) A 4 , upto A 2 = 49 and then extrapolated to A 2 = 
81. 

There are two notable features: (i) There is clearly a 
non zero coefficient at\(ri) so the equivalent of the WL 
corrections do not control the leading correction to ob in 
three dimension. The coefficients ol\ and ai are shown in 
the inset in Fig. 4. (ii) The 'low A' fit seems to work rea- 
sonably for A < 8, in the sense that p/pb ^5 2. This qual- 
itative correspondence with the Boltzmann result, even in 
the regime ao/l ^ 1, has been noticed recently 34 . 

The first issue has been explored in detail 37 by Bclitz 
and Kirkpatrick who argue that the standard WL pro- 
cesses do not exhaust the leading corrections to ob in 
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three dimension. According to them, the perturbative 
expansion for a, in a continuum model, has the form 

cr ~ a B {\ - a(M) -1 - b(k F l)- 2 log(k F l) + 0{{k F iy 2 )} 

where a and b are numerical coefficients 0(1). The WL 
argument would put a = 0, 6 = 0. 

This form for the correction beyond Boltzmann trans- 
port has apparently been observed for electron mobility 
in dense neutral gases. The detailed coefficients in this 
expression would change in a tight binding model, but 
the key result about k F l dependence should survive. 



2. Intermediate coupling 



The Boltzmann result alongwith the perturabtive 
quantum corrections is reasonable probably upto A/W ~ 
0.2 — 0.3, although numerically the fit, in the last section, 
seems to describe the resistivity even upto A/W <~ 0.75. 
The scaling regime, where localisation effects become vis- 
ible, occurs within about 10% of A c . 

Despite the correspondence of our numerical results 
with an extrapolation of weak coupling theory, there is 
no analytic framework for calculating the resistivity when 
the "small parameter" (fc^/) -1 ~ ao/l becomes 0(1). 

The paramtcr a /l is 0(1) for A/W ~ 0.5 but the de- 
viation from the Boltzmann result (at the band center) 
is only about 25%. This has been pointed out recently 
by Nikolic and Allen 34 and probably arises from a fortu- 
itious cancellation of higher order corrections. The self- 
consistent theory (SCT) of Vollhardt and Wolfle 19 serves 
as an interpolating approximation in this regime. 
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FIG. 5. The optical conductivity, at band center, for differ- 
ent L and A. Panel (a), corresponds to the moderate disorder 
regime, with a Drude form for cr(w), while (6). is for a system 
on the verge of localisation (vanishing (Jdc)- 




FIG. 6. The approach to the d.c conductivity, with increas- 
ing L in the Lt x Lt x L geometry. Lt = 6 and the disorder 
is increased from the perturbative end towards localisation. 
The chemical potential /i = 0, so n = 1. <j(L), defined earlier 
in the text, is the average of <r(w) over the interval [0, Aw], 
with Aw = 1.1/L, at Lt = 6. Inset: variation of a(L — 24) 
with A, to illustrate the rapid fall in the reference conductiv- 
ity with increasing disorder. 

Within the SCT also, the conductivity at band center 
remains within 20% of <tb for A < 8. The deviation from 
the Boltzmann result grows as we move from the band 
center to the band edge as evident in Fig. 4. 



3. Scaling regime 

The scaling regime occurs close to critical disorder, 
within about 10% of A c . The conductivity in this regime 
varies as A c — A. This regime has been extensively stud- 



ied to clarify the critical properties (see, 



e-g, 



and ref- 



erences therein). We have not used a dense enough sam- 
pling in A for discussing the critical behaviour, and our 
system sizes too are not large enough for high accuracy 
calculation of the conductivity in this regime. However, 
based on results at A = 16 and A = 17 we can bracket 
the critical point, as shown in Fig. 3. 



C. Optical conductivity 

The optical conductivity <j(uj) is of intrinsic interest 38 
and also plays a role in our method of determining the d.c 
conductivity. There are some exact results known on the 
form of the low frequency a {to) in the Anderson model. 

(i) At weak disorder, when Boltzmann transport holds, 
the optical conductivity has the Drude form, a(ui) ~ 
er(0)/(l + ui 2 t 2 ), where r _1 oc A 2 as we already know. 
For lot -c 1 this would give us a(ui) ~ er(0)(l — uj 2 t 2 ). 
(ii) When the quantum corrections to the d.c conduc- 
tivity become important the frequency dependence also 
picks up a non Drude form. In the intermediate disorder 
regime, one expects a(u>) ~ cr(0) + 0(Ty/oj/T), where 
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cr(0) already incorporates corrections beyond the Boltz- 
mann result. In this regime the conductivity rises with 
increasing frequency, for frequencies ui <C T. (Hi) At 
the critical point, where the zero frequency conductiv- 
ity vanishes, a(uj) <~ a; 1 / 3 , and in the localised regime 
(t{uo) <~ uo 2 . 

These results originally obtained through different 
techniques can be obtained in a unified way via the self- 
consistent theory of Vollhardt and Wolfle. 

Fig. 5 demonstrates the changing character of <j(uj), at 
n = 1, as we move from the Boltzmann regime (A = 
2), to strong disorder (A = 16). We show the data for 
three system sizes at each A to illustrate the explicit L 
dependence in cr(w,L). This is important for analysing 
the extrapolation needed for cr^c- 

There are two effects of changing system size: (z) the 
a{uo) profile itself can change with evolving system size, 
over the frequency range of interest, and (ii) larger sys- 
tem size allows access to (more dependable) low fre- 
quency data. 

Fig. 5. (a), the weak disorder case, reveals that the cr(w) 
profile changes perceptibly with increasing L, the changes 
being 0(5%). This implies that in our attempt to access 
d.c conductivity, the contribution arises not only from 
lowering the frequency cutoff but also moderate changes 
in the cr(u) profile. At strong disorder, Fig.5.(b), the 
profile itself does not change significantly with L and the 
key change in the ct<j c estimate comes from our ability to 
access lower frequencies. 



low frequency average a av (L), at some size L, as the 
bulk d.c conductivity? Fig. 6 illuatrates the extrapolation 
based on the sequence {L : 24,32,40,48,56,64}, at /i = 
0, moving from weak to strong disorder. 

In the weak disorder regime, the optical conductiv- 
ity is 'flat' for ui <C t -1 so if low enough frequencies 
can be accessed (given the finite size gaps), the d.c con- 
ductivity can be reasonably approximated. This is the 
feature observed at A = 2 in Fig. 6. However, in the 
WL region and beyond, ct(uj) has non trivial frequency 
dependence at low u, as evident in Fig.5.(b). The cor- 
responding low frequency average has significant L de- 
pendence. Since <t(uj) ~ er(0) + 0(y/uj), the low fre- 
quency average a av (L) ~ cr(oo) + 0(1/VZ). The data 
at A = 12,16 show a reasonable fit to the square root 
form. The much stronger frequency dependence in the 
strong disorder regime makes a size dependent study im- 
perative. We provide a discussion of the extrapolation 
scheme in an appendix. 

These results illustrate the work involved in access- 
ing the d.c conductivity, particularly in the regime of 
strong disorder, where a small L calculation (at L = 16 
say) might overestimate the conductivity by a factor of 4. 
This discrepancy worsens as A — > A c and a systematic 
study of size dependence is vitally important. 



E. Effect of the transverse dimension 



D. Large L extrapolation 

How important is the large L extrapolation to access 
the d.c conductivity, i. e, what is the error if we treat the 
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FIG. 7. Effect of transverse dimension on the large L ex- 
trapolation. Carrier density n = 1. The conductivity ob- 
tained by extrapolating the 6 x 6 x L results are shown as 
circles on the y axes. 



All the results quoted till now have been obtained via 
extrapolation on a 6 x 6 x L geometry. The 6 2 cross 
section was chosen to allow large L to be accessed. How- 
ever, it is important to quantify the error involved in 
chosing a specific transverse dimension Lt- To this end 
we studied the low frequency average a av in a sequence 
Lt x Lt x L with L T = 2, 4, 6, 8, 10, 12 and L = 10, 20, 30, 
for A = 2, 10 and 16, and n = 1. The averaging interval 
Auj was scaled as 1/(L 2 Lt) in all geometries. 

Fig. 7 shows a av with respect to for the sequence 
Lt specified above. Panel (a), shows the weak disorder, 
A = 2 result. Beyond Lt — 4 all the curves seem to 
converge to a ~ 0.70 for L — ► oo. The extrapolation 
from Lt = 6, obtained using L upto 64, is shown as a 
circle on the y axis, and is ~ 0.68. 

For the strong disorder case, Fig.7.(b), the extrapola- 
tion for Lt = 6 is shown to be ~ 0.001, while the large 
Lt data, using L upto 30, suggests that the asymptotic 
value could be larger, <~ 0.002. This suggests that 'small' 
Lt somewhat underestimates the conductivity (remem- 
ber Lt — 1 is one dimensional, so completely localised), 
while finite L overestimates the conductivity. Except 
very close to the MIT these errors are small for the sizes 
we use and, as verified by the phase diagram, even the 
critical point is located to within 5%. 
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IV. SCATTERING FROM MAGNETIC 
DISORDER 

A. Global features 

The effect of weak magnetic scattering on transport 
is quite similar to that of potential scattering. The ef- 
fect is contained in the Born scattering rate, r" 1 oc 
N(ep)J' 2 S 2 , and the weak coupling resistivity p(J',n) 
varies as ~ &i (n) J' 2 , where b\(n) is a density dependent 
coefficient. However, even at moderate coupling, J 1 ~ 2, 
new effects begin to show up in cr(/i). The conductivity 
at half-filling, n ~ 1, gets suppressed more quickly than 
would be guessed based on the Born argument. This de- 
viation, and its evolution with increasing J', arises from a 
fundamental difference between potential scattering and 
magnetic scattering on a 'Kondo lattice'. 

There are in fact two main differences that show up be- 
yond weak coupling. These are visible when we compare 
Figs. 8-10, with Figs. 1-4. (i) The conductivity in the po- 
tential scattering case decreases monotonically (at fixed 
n) with increase in disorder, vanishing at A c (n), while in 
the magnetic scattering case, at a generic density, the 
conductivity is finite even as J' — > 00. The resistiv- 
ity 'saturates' and there is no metal-insulator transition 
with increasing J', except in a narrow density window. 
(ii) The band center, n ~ 1, is of no particular signifi- 
cance in the Anderson problem, except A c being largest. 
In the J' problem the response for n ~ 1 is dramatically 
different from that in the rest of the band. There is an 
MIT at J' w 5. These differences can be understood from 
an analysis of the strong coupling end. 

For J 1 It ^> 1 it is useful to choose a local quantisation 
axis at each site, for the electrons, parallel to the orien- 
tation of the spin Sj. The coupling J'Sj acts as a strong 
local Zeeman field on the electron. Suppose the hopping 
term were absent. The two local eigenfunctions at each 
site would have spin projections parallel and antiparallel 
to J'Si, with energy /2 respectively. The zero hop- 
ping problem leads to N fold degenerate levels at ± J'/2. 
The 'gap' J' plays a key role at strong coupling. The 
presence of hopping generates a degenerate perturbation 
on the locally aligned states (say), and the electrons can 
now 'hop' with an amplitude that depends on the orien- 
tation of nearest neighbour spins. This mechanism has 
been extensively discussed in the context of the double 
exchange model. 

The mixing introduced by 'hopping' broadens the two 
levels into bands. For J' > W/2 the broadening due to t 
does not fill the gap, and the system is insulating at n = 
1. For J' below this critical value, J c , say, the DOS at 
band center is suppressed but finite, and the resistivity is 
still very large. In summary, the strong coupling physics 
of (incipient) band splitting controls the resistivity close 
to band center, and creates an essential difference, in 
terms of J' and n, with respect to standard Anderson 
localisation. 
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FIG. 8. (a). Conductivity as a function of Fermi energy, 
(6). density of states, for different values of J', in the case of 
pure magnetic scattering. 
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FIG. 9. Dependence of conductivity on carrier density for 
varying J'. The conductivity scale is logarithmic. 
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FIG. 10. Resistivity variation with J', for several electron 
density, from the perturbative limit to double exchange 
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The saturation in p(J') with increasing J', over most 
of the band, occurs because the effect of large J' is ab- 
sorbed mainly in the band splitting. The effective disor- 
der seen by the electrons comes from fluctuations in the 
hopping amplitudes, explained in the next section, and 
these are 0(t). The ratio of fluctuation to mean hopping 
is moderate, so the large J' limit leads to a 'dirty metal' 
but no metal-insulator transition. This is unlike the An- 
derson problem where the electrons scatter off potential 
fluctuations whose amplitude grows with increasing A. 

From the data in Figs. 8-10 we can now identify the 
different transport regimes. 

B. Transport regimes 

1. Weak coupling: J' /W <§C 1 



exception is the vicinity of n = 1, and the lower edge 
of the band. This suggests that the correction to the 
Born resistivity is negative for n away from n = 1, and 
changes sign as n — > 1. Fitting the data to p(J',n) = 
bi(n)J' 2 + 62(71) J' 4 , the coefficient b2(n) illustrates the 
crossover from saturation to escalation, as we move across 
the band. Fig. 11 shows the fit to this form and the coef- 
ficients are shown in the inset. The 'Born' coefficient is 
positive throughout the band, without significant density 
dependence in the density interval shown. The quartic 
coefficient changes sign, from positive to negative, as n 
is lowered from 1.0 to 0.74. 

A confirmation of saturation or escalation cannot of 
course be obtained from a low order expansion in J' 2 , 
but even the 'perturbativc' coefficient provides a hint of 
strong coupling physics. It also suggests a smooth evolu- 
tion from weak to strong coupling. 



The magnetic scattering rate T s is proportional to 
N(€f)J' 2 , and the weak coupling resistivity should be ex- 
pandable in r s . The lowest order term is well known, cor- 
responding to Born scattering, with p(J',n) <~ bi(n)J' 2 . 
The density dependence is similar to that for potential 
scattering. Assuming p(J') to be analytic in r s , i.e ig- 
noring possible log corrections etc, Fig. 11 shows a fit of 
the form p(J', n) ~ bi{n)J' 2 + b 2 (n)J' 4 to the low J' re- 
sistivity. The J' 2 character dominates upto J' ~ 2, as 
one can see also in the a(n) plot in Fig. 9, beyond which 
the quartic term becomes important. We do not know if 
the coefficient of the quartic term has been analytically 
calculated, but the sign of this term is crucial, and is 
density dependent, as we discuss next. 



2. Intermediate coupling: J' /W ~ 0(1) 

As is obvious from the data in Fig. 9-10, the resistivity 
saturates with increasing J', over most of the band. The 



20 






FIG. 11. Fit to the weak coupling resistivity of the form 
p(J',n) ~ b\{n)J' 2 + &2(n)J' 4 . The symbols are actual data 
and the firm lines are fits. The inset shows the n dependence 
of the coefficients bi and &2- Note the sign change in &2- 



3. Double exchange: J' /W — ► 00 

In the double exchange limit the J' scale acts as a 'con- 
straint' on the electron spin orientation and no longer 
directly affects physical properties, the only effect is to 
renormalisc the chemical potential. The mapping of the 
J'/t — > 00 problem to a 'spinless fermion' problem with 
hopping dependent on nearest neighbour spin orienta- 
tion has been widely discussed 39 . Transformed to spin- 
less fermions, which correspond to original electron states 
with spin projection 'locked' parellel to the local quanti- 
sation axis, Sj, the Hamiltonian becomes: 



H 



(ij) 



(ij) 



(ij) 



(5) 



The tij being the spin orientation dependent hopping am- 
plitude specified earlier. We can split it into the mean 
(uniform) hopping amplitude, i, and the fluctuation SUj. 

In the 'extreme' paramagnetic phase of this model, the 
distribution of hopping integrals is exactly known. The 
spins are independently distributed on a sphere so the 
can be worked out. There is no obvious small param- 
eter, since both the mean value of hopping, f, as well 

as the fluctuation, At = (St 2 ), are cx t. However, 

the ratio At/i w 1/3. Numerical work by Li et al. 2A 
had demonstrated that less than 0.3% of states in the 
band are localised under this condition. It was not clear 
whether the resistivity at the band center, n = 0.5, could 
be described within a Boltzmann approach. Narimanov 
and Varma 40 have demonstrated that the mean free path 
emerging from the Boltzmann calculation is l/ao > 8 so 
the method is self-consistent. 

It seems now that despite the localisation effects as 
n — ► 1 and n — > 0, resistivity over much of the band 
can be understood within a effective 'weak coupling' ap- 
proach. The resistivity is w (0.1 — 0.2)pMott at the band 
center according to our calculation. The resistivity is also 



9 



'particle-hole' symmetric, now within the lower band, but 
notice that this is cleanly visible only at very large J' ft. 

4- Virtual orbital mixing: large finite J' 

As we move to finite J' from the double exchange limit, 
the two bands still remain split (down to J'/t s=s 5) but 
there is a virtual admixture that is introduced. To access 
properties in this regime we need to use a two orbital 
formulation, with the orbital energies still separated by 
a large gap ~ J'. The chemical potential remains in the 
lower band. The two orbital model, written in terms 
of electronic states with local quantisation axis, has the 
form: 

ij i i 

We have not seen a Boltzmann calculation of transport 
in this regime, but using the two orbital formulation it 
might be possible to set up such a scheme. The resistivity 
decreases as we move down from large J', so using the 
correct 'basis' the transport may be accessible within a 
Boltzmann approach (since the double exchange limit is 
itself so accessible). The increase in conductivity, da(J'), 
as we move to lower J', is found to be proportional to 
1/J'. A perturbative correction to the large J' result, 
within a diagrammatic scheme yields the same answer. 



We do not know if any analytic approaches have been 
explored in this localisation problem. This regime would 
be relevant to the low doping magnetic semiconductors, 
where there is also the possibility of carriers trapping into 
spin polaronic states. 



C. Optical conductivity 

The optical conductivity confirms the trends seen in 
the d.c conductivity. Fig. 12 shows <j{uj) at J' = 1 and 
J' = 4, weak and 'strong' coupling respectively, at the 
center (n = 0.5) and edge (n = 1) of the lower band. 

At weak coupling, over the frequency range shown, 
o~(u>) is larger at n = 1, the center of the full band, com- 
pared ton = 0.5. The scattering rate is T s = 2nN(eF)J' 2 
which at the band center is ss 0.75. If a(uj) follows the 
Drude form then <7(r s )/cr(0) should be ~ 0.5, which is 
consistent with Fig. 12. (a). By J' = 4, the trend has 
reversed. The n = 1 case is almost insulating, with 
er(0) — > 0, while the conductivity at the (lower) band cen- 
ter is finite and essentially flat on the scale considered. 
This trend gets amplified as we go to even larger J'. 



V. COMBINED STRUCTURAL AND MAGNETIC 
DISORDER 

A. Global features 



5. Behaviour near band tails 

Spin disorder by itself cannot localise states in the cen- 
ter of the band, since At/t is not large enough. However 
at the band tails, i.e, n close to 1 or 0, the kinetic en- 
ergy is small and a small fraction of states can still be 
localised. As we have indicated earlier, this is < 0.3% 24 
of the total number of states for J' — > 00. 




Q) to 



FIG. 12. Optical conductivity: (a). Drude response at both 
n = 0.5 and n = 1.0 at weak coupling, and (b). strong scat- 
tering at n = 0.5 and alomst 'insulating' response at n = 1.0 
at strong coupling. 



In the presence of both structural and magnetic dis- 
order it is not possible to show the full density depen- 
dence of transport properties compactly, so we provide 
two generic 'cross sections' in Fig. 13 at n = 0.26 and 
n = 1.00. In addition to the effects already noted for 
potential scattering and magnetic scattering, there are 
several novel features that arise. 

(a) . For weak J' and moderately large A, magnetic 
scattering weakens localisation effects, as evident from 
the intermediate A small J' data in Fig. 13. 

(b) . At even larger A, where the system would have 
been Anderson localised, magnetic scattering converts the 
insulator to a metal. The critical disorder A c (n) shifts 
to a larger value A c (n, J'), see phase diagrams in Fig. 14. 

(c) . In contrast to purely magnetic scattering, where 
the resistivity typically 'saturates' with increasing J', in 
the presence of structural disorder the system can go in- 
sulating with increasing J'. 

(d) . The 'additivity' of magnetic and structural scat- 
tering holds only over a very limited range in A and J', 
Mathiessens rule generally does not hold. 

The major features, above, can be easily motivated 
after we write down the different effective models of scat- 
tering in the various transport regimes in the problem. 
Some of this has been discussed earlier by us 27 , so we will 
discuss mainly those aspects of the problem which have 
not been covered earlier. 
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B. Transport regimes 

The parameter space of the problem is large, involving 
n — A — J', and it is convenient to first identify distinct 
density ranges and then classify the transport/scattering 
mechanisms. The roughly distinct density regimes are 
the following. 

(1) . The wide 'mid band' region 0.1 < n < 0.9, of which 
n = 0.26 in Fig. 13. (a) is typical, and we have discussed 
the n = 0.5 case earlier 27 . 

(2) . n — > 1, where the response is similar to the mid 
band region at weak J', but the large J' response is dis- 
tinctive, Fig.l3.(6). 

(3) . n — > 0, where structural disorder and magnetic 
scattering readily leads to localisation. 



1. Generic density: the mid band region 

Let us consider this typical density regime first. For 
generic densities, 0.1 < n < 0.9, say, there are tentatively 
five different transport regimes in the problem. These 
are: (i) Both A and J 1 small: the weak scattering regime, 
where the effect of structural disorder and magnetic scat- 
tering are perturbative and additive, (ii) Moderate A 
and small J': spin flip correction to weak localisation. 
The A dependence shows WL corrections and spin flip 
scattering weakens the WL correction. (Hi) Large A, 
~ A c , and small J': spin dephasing driven insulator- 
metal transition (IMT). (iv) J' /t — > oo, with varying A: 
the disordered double exchange (DE) limit, (v) J'/t 3> 1 
but finite, and moderate to large A: the intermediate 
coupling 'metal'. 



10 20 30 




FIG. 13. Global behaviour of the resistivity with vary- 
ing structural disorder (A) and magnetic coupling (J'): 
(a), n — 0.26 typical of most of band, and (b). n = 1.0 which, 
at large J', corresponds to the upper edge of the lower band 
and has its own distinct transport response. 
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FIG. 14. Putting together different constant density 
cross-sections to create a 'global' insulator-metal phase di- 
agram for the A + J' problem. The densities are marked in 
the panels. M stands for a metallic phase while I is insulating. 
The bounding curves can be viewed as A c (n, J'). Notice the 
log scale on J'. 



(i) When A and J' are both small the transport can 
be understood in terms of additive Born scattering, with 
the net scattering rate, T(A, J') a\ A 2 + bi J' 2 , and the 
resistivity p oc T(A, J'). The 'window' describing this 
regime is roughly J' < 3 and A < 4. The resistivity 
in this regime is p < O.lpMott, *-e, below 100 /if2cm, 
say. This corresponds to the bottom left hand corner in 
Fig. 13. (a) and, as our earlier data showed 27 , Mathiessens 
rule holds. 

(ii) At larger A remaining at small J', as the WL cor- 
rections show up, spin flip scattering 22 ' 23 of the electrons 
by the random magnetic moments reduces the localising 
effect of structural disorder, i.e, dp/dJ'\ 7lt \ < 0. Just as 
inelastic scattering weakens quantum interference by in- 
troducing decoherence, spin flip scattering leads to spin 
decoherence. We have quantified the A and J' depen- 
dence of the effect in the earlier paper 27 . 

(Hi) At even larger disorder, A > A c , where the J' = 
system would have been Anderson localised, spin flip 
scattering opens up a metallic window. The structural 
disorder needed for localisation shifts to a larger value, 
i.e, dA c /dJ'\ n > 0. This effect in visible in all the panels 
in Fig. 14. 

(iv) Now consider the DE limit, J' — ► oo. As we 
have discussed in Section IV, the form of the resistivity 
p( J', A = 0) arising from 'magnetic disorder' at large J' 
is very different from what one observes in p(J' = 0, A) 
at large A. This is because J' contributes to both 'band 
splitting' and effective disorder, and the effective disor- 
der saturates as J'/W — > oo with J' controlling only the 
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band splitting. The presence of structural disorder in the 
J' — > oo problem strongly enhances the resistivity and lo- 
calising tendency. Using the transformations used in the 
previous section: 

ij i 

= toihi + E 7 ^ + E e ^ 7i (7) 

ij i 

The localisation properties of this model have been stud- 
ied by Li et aZ. 24 , although they did not calculate the 
resistivity. The 'hopping disorder' by itself localises less 
than 0.5% of the states in the band. On adding structural 
disorder the mobility edge moves inward, with localisa- 
tion of the full band occuring at A/t ~ 11.5, which can 
be approximately understood from the roughly 30% band 
narrowing due to spin disorder. 

(v) Finally the regime with large but finite J' and 
strong structural disorder. We have seen that the effect 
of small J' at strong structural disorder can be qualita- 
tively understood in terms of the spin dephasing effect on 
Anderson localisation. However, the small J' behaviour 
with dp/dJ' < quickly leads to a minimum and then a 
regime with dp/dJ' > 0. Such behaviour can be viewed 
as an extension of the J' 2 term seen at weak disorder, but 
it is more fruitful to approach the effect from the strong 
coupling DE end, as we do below. Transforming to the 
usual local spin quantisation frame and retaining both the 
parallel and anti-parallel electron states, we have: 

ij i i 

The major source of disorder is still e i; with additional 
contribution from the 6t"f. The orbital mixing effect 
of 'off diagonal' couplings, either in terms of mean am- 
plitude or fluctuations, is regulated by the large en- 
ergy denominator J' . Although the 'reference' problem, 
J' — > oo, is not analytically tractable in the presence of 
structural disorder, it can be shown that orbital mixing 
generates a correction to conductivity <~ 0(1/ J'). 

2. Half -filling: n — » 1 

For n — > 1, the effects at small J' are similar to (i) — 
(Hi) at generic densities, discussed above. This is borne 
out by the behaviour of p(J',A) in Fig. 13. (b) and the 
phase diagram in Fig. 14. (d). At large J', however, the 
system always goes insulating, see Fig. 13. (b), as we have 
discussed in Section IV as well. This effect is obviously 
due to the band splitting induced by large J' and the 
vanishing DOS, N(ep), at n = 1.0. Thus, for n = 1.0, the 
metallic phase is bounded both in A and J', Fig. 14. (d). 

There is however an interesting and possibly unex- 
pected feature in Fig. 14. (d) for J > 5.0, where the A = 
system becomes insulating. We may have imagined that 



introducing structural disorder in this system would en- 
hance localisation. This however is not true, and struc- 
tural disorder actually 'metallises' the reference band 
split state, and the critical J' needed for localisation in- 
creases in the presence of structural disorder. 

The origin of the effect above lies in the 'band broad- 
ening' effect of structural disorder. The A = problem 
had a narrow (vanishingly small) gap in the DOS, and 
the presence of structural disorder creates finite DOS at 
the Fermi level, effectively closing the gap. Since the net 
disorder arising from the random spins and the structural 
disorder is relatively weak the finite DOS seems sufficient 
to lead to a metallic, albeit highly resistive, phase. For 
A ~ 4 and J' <~ 5, the resistivity is roughly Q.bpMott- 
As A becomes large, or J' becomes large, this metal- 
lic window is lost due to the effects either of A driven 
localisation or band splitting. 

3. Very low density: n — ► 

The case of n — > 0, for example n = 0.01, say, is un- 
fortunately hard to access with control for the system 
sizes that we have used. We expect that the small J' 
behaviour will be similar to that in the rest of the band, 
with enhanced resistivity (due to the low carrier density) 
while the behaviour for J' — > oo will be similar to that 
for n — > 1 (due to the particle-hole symmetry within the 
lower band, < n < 1). Localisation in the n — > limit, 
we believe, is better explored via transfer matrix methods 
due to the large accessible size. 

VI. CONCLUDING REMARKS 

In this paper we have presented controlled results on 
electron transport in the background of arbitrary struc- 
tural and spin disorder and provided a framework within 
which the data can be analysed. We benchmarked our 
Kubo formula based method in the standard problem of 
potential scattering and Anderson localisation. We then 
explored the distinct transport regimes that arise in the 
case of pure magnetic scattering, as well as the combined 
effect of structural and magnetic disorder. In contrast to 
the effect of only structural disorder (where the resistiv- 
ity 'escalates' with increasing disorder) or only magnetic 
scattering (where it 'saturates' with increasing disorder) 
their combined action can lead to non monotonic depen- 
dence and novel transport regimes. The method devel- 
oped in this paper can be directly taken over in calculat- 
ing the resistivity in the presence of annealed disorder, 
where accessible system sizes rarely exceed ~ 10 3 , and 
has been extensively used by us in Monte Carlo studies 
of several "disordered" electron systems. 

We acknowledge use of the Beowulf cluster at HRI. 
S.K. gratefully acknowledges support by the Deutsche 
Forschungsgcmcinschaft through SFB484. 
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VII. APPENDIX: EXTRAPOLATION FOR THE 
D.C CONDUCTIVITY 

It is known that at sufficiently low frequency the op- 
tical conductivity in the 3D Anderson model follows a 
simple power law 19 with the exponent depending on the 
strength of disorder and electron density. This means 
that we can write: cr(u)) = A + Buo a , where the co- 
efficient A > is <Jdc, while the next term gives the 
leading low frequency correction. This form captures 
all the broad regimes in 3D. For example, at half-filling: 
(i) low disorder, Born scattering, A <C A c , gives B < 0, 
a = 2, (ii) moderate disorder, weak localisation correc- 
tions: B > 0, a = 1/2, (Hi) critical disorder, A = A c : 
A = 0, B > 0, a — 1/3, and (iv) localised phase, A > A c : 
A = 0, B > 0, a = 2. 

The form for <j(u>) fixes the form for <r av (Aui, L). Set- 
ting Auj ~ 1/L, and using the form for a(oj) above, we 
obtain the three parameter form for a av (L): a av (L) ~ 
A + r^j\ L~ a . The extrapolation is a least square three 
parameter fit to our L dependent data, and has enough 
flexibility to cover all disorder regimes. 

In practice, a simpler two parameter extrapolation also 
works reasonably as long as one is in the metallic phase, 
even close to the metal- insulator transition: <j av (L) <~ 
A + ^L- 1 ' 2 This derives from a(w) ~ A + Buo 1 ' 2 . 
Using this restricted fitting function, the MIT can be 
roughly located when A, the d.c conductivity, falls below 
a preset limit, 10~ 6 say, (which in absolute units is a 
resistivity ~ 10 4 puott)- Having located the transition 
approximately, the more elaborate three parameter fit 
can be used to confirm the metallic/insulating character 
on two sides of the critical point. Our MI phase diagrams 
are constructed using this strategy. 
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